Large thermoelectric figure of merit for 3D topological Anderson insulators via line 

dislocation engineering 
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We study the thermoelectric properties of three-dimensional topological Anderson insulators with 
line dislocations. We show that at high densities of dislocations the thermoelectric figure of merit ZT 
can be dominated by one-dimensional topologically-protected conducting states channeled through 
the lattice screw dislocations in the topological insulator materials with a non-zero time-reversal- 
invariant momentum such as Bii_xSbx. When the chemical potential does not exceed much the 
mobility edge the ZT at room temperatures can reach large values, much higher than unity for 
reasonable parameters, hence making this system a strong candidate for applications in heat man- 
agement of nano-devices. 

PACS numbers: 73.50.Lw, 71.90.-(-q, 72.20.Pa 
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Introduction. The recent crisis of heat management in 
nano-devices, which has lead to a lack of progression in 
clock speeds of charge-based logic devices, has intensified 
the interest in efficient thermoelectric materials. In the 
last decade there has been a lot of research both the- 
oretical [IHS] and experimental [6l-(8] to create efficient 
thermoelectric nano-devices. The efficiency of such ma- 
terials, which must be both p-type and n-type, is de- 
termined by a balance to convert charge flow into effi- 
cient heat transport as well as maintaining a temperature 
gradient between the device and the heat sink. Among 
the most well known thermoelectric materials in present 
day commercial applications one finds Bi2Te3, PbTe, and 
PbSb. These type of insulators or semiconductors have 
been identified recently as topological insulators (TI) [51 
which exhibit protected delocalized surface states. 

In the two-dimensional version of the TIs, the quan- 
tum spin Hall systems [101 III]j these protected edge 
states contribute to the thermoelectric efficiency but do 
not enhance it dramatically beyond present day materi- 
als whose efficiency parameter, ZT (see below), is of 1 or 
slightly less [2.;. On the other hand, for three-dimensional 
(3D) topological insulators with a non-zero time-reversal- 
invariant momentum [H] (TRIM) it has been shown the- 
oretically that one-dimensional (ID) topologically pro- 
tected modes can exist in the bulk propagating through 
certain line dislocations [T3]. These type of ID quan- 
tum modes have been recently attributed to the results 
of recent experiments on Bi2Se3 [l4) , fT5 ] . 

Here we explore the idea of using these ID topologi- 
cally protected modes to significantly increase the ther- 
moelectric efficiency of materials such that Bii_xSbx, see 
Fig. [T] The basic premise of the proposal is to introduce, 
through growth engineering, a finite density of screw dis- 
locations. This would induce disorder in the bulk leading 
to a reduction of the thermal conductivity, Anderson lo- 
calization of bulk states, and an increase of the conductiv- 




FIG. 1: (Color online) A sketch of 3D topological insulator 
with lattice dislocations propagating through it. These dislo- 
cations are topologically-protected ID modes which are per- 
fectly conducting. The inset shows the details of the screw- 
type dislocation considered. 



ity and thermopower contributions from these ID states. 
This combination of factors, as shown below, leads to 
a dramatic enhancement of the figure of merit efficiency 
for thermoelectrics, ZT, beyond its present value for bulk 
materials. For reasonable parameters we estimate ZT to 
reach ^ 6 at room temperature. 

Figure of merit. The performance of thermoelectric 
devices is determined by the dimensionless figure of merit 
ZT defined as 



ZT = 



aS^T 



(1) 



where a, S, k, and T are the electrical conductivity, 
Seebeck coefficient (or thermopower) , thermal conductiv- 
ity, and absolute temperature, respectively. The thermal 
conductivity generally has two contributions: electronic 
and phononic, k = Kg^-'rKpii. A large effort within thermo- 



electric studies is devoted to finding materials with high 
ZT near room temperature. The challenge is that effects 
that increase a o\ S usually accompanies an increase on 
K and vice versa. In particular, an important challenge is 
to reduce the phononic thermal contribution to n which 
can limit severely ZT . Doing so without reducing at the 
same time the aS'^T factor is a key challenge. 

Our proposal is to exploit the ID topological protected 
channels formed at the dislocations of a TI with non-zero 
TRIM [13] to circumvent this challenge. As more dislo- 
cations are introduced in the system the mean free path 
of the phonons is reduced while at the same time the con- 
ductivity contribution from the ID channels is increased 
as mentioned above. By requiring that the bulk contri- 
bution is also reduced by Anderson localization of the 
bulk states, through the dislocations themselves or other 
co-doping, the ZT factor will be dominated exclusively 



by the ID channels and a large ZT value can be obtained 
|16j . Below we present our theoretical estimate for this 
high ZT using reasonable estimates of the different pa- 
rameters without seeking the best possible scenario but 
estimating instead a reasonable expectation of the pro- 
posed system. 

Within linear response theory [T7j the electric {j^) 
and thermal (j'^) currents are given by linear combina- 
tions of the chemical potential and temperature gradi- 
ents: f/e = LoV^i -I- Li(Vr)/r and j« = -LiV/i - 
L2(Vr)/T, where e is the electron charge. From these 
equations, using Onsager relations, one can find that 
the electrical conductivity a = e^Lg, thermopower S = 
— Li/(eTLo), and electronic thermal conductivity Kg = 
{LqL2 — L\)/{TLq). Then according to Eq. ^ in terms 
of transport coefficients La the figure of merit takes the 
form [18j : 



ZT = 



{Ll 



snL, 



^i^)(L| + snLl^) - {L\ + snLY'f + KpH{Ll + snLlD)T' 
I 



(2) 



where n is the density of topologically-protected lat- 
tice dislocations [13], s is the cross-sectional area of 
the device transverse to the transport direction, and 
Kp/i = -^Vphlph is the phonon contribution to the ther- 
mal conductivity. This expression for Kph is applicable 
at room temperatures. In Eq. (pi) it is assumed that the 
transport coefficients have bulk and ID channel contri- 
butions La = L^a + snLa^ , whcrc sn gives the number of 
the topologically protected lattice dislocations conduct- 
ing perpendicular to area s. 

Next we take Iph in the limit of high density of ran- 
domly located dislocations, such that for high enough dis- 
location density n, Iph is diminished by phonon scattering 
from these dislocations. The phonon specific heat Cph at 
room temperature can be estimated as Sriikg, where the 



number of ions per unit volume is 



4 • 10 



29, 



In (Bii_2;Sb:r)2Te3 and Bi2(Tei_j^Sej,)3 compounds the 
phonon velocity Vph = 1500 m/s and Debye tempera- 
ture 9u = 142 K [19] . We also take into account in 
our calculations the fact that the bulk becomes an amor- 
phous media at very high n due to disorder which leads 
to the saturation of Kph- Therefore, Kph ranges from 
lWm~^K~^ for pure bulk with no dislocations |20j to 
Kph « O.OlWm^^K^^ at average distances between dis- 
locations d ^ 3 nm. It is clear from Eq. ^ that at large 
densities of lattice dislocations, n, the contribution to ZT 
mostly comes from ID channels and in the limit equals 
to ZT of one perfectly conducting ID wire: 



lim ZT 



iL\ 
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(Ln 



(3) 



To make estimates of the relative contributions of the 
ID channels and bulk we model the topological Anderson 
insulator system as a semiconductor with one valence and 
one conduction band and one ID-state corresponding to 
a perfectly conducting lattice dislocation (in general, one 
ID-state per each dislocation), see Fig. [2] (a). We take the 
gap between conduction and valence bands in the bulk to 
be A (—A < E < 0) and the transmission coefficient in 
the ID channel is assumed to be unity T{E) = 1 for sim- 
plicity Generally A = 0.15 eV for Bi2Te3 and A = 0.3 
eV for Bi2Se3 we use the former for our estimates. The 
bulk states near the edges of the band are localized with 
a mobility edge at E^ ', this Anderson bulk localization is 
due to the high density of dislocations or can be induced 
by doping with non-magnetic impurities. 

Employing the Landauer-Buttiker formalism [2ll [22] 
we can write the expressions for the transport coefficients 
in the ID channel: 



ID 



I 



r{E)f{E){E~^irdE, 



(4) 



where I is the length of the sample in the growth direction 
(length of ID channel) which has the upper limit of the 
inelastic coherence length {lin ^ 1/^rn), his a Planck con- 
stant, and f'{E) = df/dE with / = l/(e(-s-^)/('=«^) + l) 
being the Fermi function. Here the integration over ener- 
gies extends from —A to while the chemical potential, 
which can be changed by an external gate, we restrict to 
be in the gap, —A < /i < 0, or within the conduction 
band, /i > 0, but not much above the mobility edge. The 



(a) E 



(b) 




15 


ZT 


/ 


5 




/ 


n 




^li/{kBT) 



Eq. ([5| gives 



3 -2 


-1 


1 2 3 


ZT 


/ 


/./(fcflT) 



! -2 -1 1 2 3 

FIG. 2: (Color online) (a) Schematic band structure of the 
bulk and ID states, (b) Figure of merit ZT for the contribu- 
tion from topologically-protected ID channels as a function of 
chemical potential ^ (in units ksT) at temperature T = 300 
K. (c) Pure bulk contribution to ZT as a function of /x. 
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FIG. 3: (Color online) The contour plot of figure of merit ZT 
as a function of chemical potential /i (in units ksT) and the 
density of lattice dislocations n at temperature T = 300 K 
for sample's length I — 1/xm. The maximum value is higher 
than 6. 



latter case is of the most interest, since in this case ZT 
is the largest, see Fig. [2](b). 

To estimate the bulk contribution to the transport co- 
efBcients Lq we assume the bands to be parabolic and 
use Boltzmann equation in the relaxation time approxi- 
mation, 



l: 



/•OO 

JE,„ 



(5) 



Here v is electron velocity, D{E) is the density of states, 
and we approximate relaxation time r to be independent 
of energy. The extended and localized states are sepa- 
rated by the mobility edge E^ which is measured from 
the bottom of the conduction band, see Fig. [2] (a), and 
for further estimates we assume Em = 0.05 eV. Then 



L' 



2^/2^ 



2/13 



:^cT"+3/2 



dx 



a;"(x-|-/i/T)3/2e^ 
(e^ + 1)2 ' 



(6) 
where c is the number of the carrier pockets and we take 
c = 1, m* — 0.02?7T,e is the effective mass, and tempera- 
ture is measured in units of ks- We estimate t = 10^^"* 
s. These L^ are then substituted into Eq. pi to give ZT 
which is shown in Fig. |3] for a range of densities n and 
chemical potential ji at room temperature. The max- 
imum, for the estimated parameters, is higher than 6, 
which makes these systems quite unique even if some- 
thing of that order can be reached. 

We note that, if the impurities are non-magnetic, the 
transition to a bulk Anderson insulator should not de- 
stroy the topologically protected ID states, because the 
time reversal invariance is not broken. Also, note that by 
increasing number of dislocations or non-magnetic impu- 
rities one can come closer to the ID wire limit for ZT 
and reach even greater values. These densities are how- 
ever more unlikely to be achieved and can lead also to 
tunneling between the channels which could lead to an 
opening of a gap and localization of the protected states. 
This is the reason for our choice of an upper limit of the 
ID-channel density of n ^ lO^^m"^ corresponding to a 
typical spacing of ^ 3 nm. 

Summary. We have studied the thermoelectric 
properties of 3D topological insulators with non-zero 
TRIM which contain many line dislocations possessing 
topologically-protected perfectly conducting ID states. 
We have shown that in principle this system can have 
very high figure of merit, ZT ^ 10, and predict that 
increasing the number of dislocations in a TI film will 
exhibit an increase in ZT. 

This work was supported by NSF under Grant No. 
DMR-0547875 and Grant No. 0757992, by the Research 
Corporation Cottrell Scholar Award, and by the Welch 
Foundation (A-1678). J. S. and S. M. thank the KITPC- 
Beijing for their kind hospitality during which part of this 
work was done (activity name "Progress in Spintronics 
and Graphene Research"). 



[1] Y. Dubi and M. Di Ventra, arXiv:0910.0425. 

[2] R. Takahashi and S. Murakami, Phys. Rev. B 81, 161302 

(2010). 
[3] T. Markussen, A. -P. Jauho, and M. Brandbyge, Phys. 

Rev. Lett. 103, 055502 (2009). 
[4] S. Mukerjee and J. E. Moore, Appl. Phys. Lett. 90, 

112107 (2007). 
[5] P. Ghaemi, R. Mong, and J. E. Moore, arXiv: 1002. 1341. 
[6] R. Venkatasubramanian, E. Siivola, T. Colpitts, and 

B. O'Quinn, Nature 413, 597 (2001). 
[7] H.-K. Lyeo, A. A. Khajetoorians, L. Shi, K. P. Pipe, 

R. J. Ram, A. Shakouri, and C. K. Shih, Science 303, 



816 (2004). 
[8] Y. Zhang, C. L. Hapenciuc, E. E. Castillo, T. Borca- 

Tasciuc, R. J. Mehta, C. Karthik, and G. Ramanath, 

Appl. Phys. Lett. 96, 062107 (2010). 
[9] Z.-L. Xi and S.-C. Zhang, Physics Today 63, 33 (2010). 
[10] B. A. Bernevig, T. L. Hughes, and S.-C. Zhang, Science 

314, 1757 (2006). 
[11] M. Konig, S. Wiedmann, C. Brune, A. Roth, H. Buh- 

mann, L. W. Molenkamp, X.-L. Qi, and S.-C. Zhang, 

Science 318, 766 (2007). 
[12] L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 

106803 (2007). 
[13] Y. Ran, Y. Zhang, and A. Vishwanath, Nat. Phys. 5, 298 

(2009). 
[14] J. G. Checkelsky, Y. S. Hor, M.-H. Liu, D.-X. Qu, R. J. 

Cava, and N. P. Ong, Phys. Rev. Lett. 103, 246601 



(2009). 
[15] Y. S. Hor, A. Richardella, P. Roushan, Y. Xia, J. G. 

Checkelsky, A. Yazdani, M. Z. Hasan, N. P. Ong, and 

R. J. Cava, Phys. Rev. B 79, 195208 (2009). 
[16] J. Sinova, unpublished. 
[17] N. W. Ashcroft and N. D. Mermin, Solid State Physics 

(Sanders College Publishing, 1976). 
[18] M. P. Marder, Condensed Matter Physics (John Wiley & 

Sons, Inc., 1976). 
[19] K. Yokota and S. Katayama, Jpn. J. of Appl. Phys. 12, 

1205 (1973). 
[20] A. Pattamatta and C. K. Madnia, Int. J. of Heat and 

Mass Transfer 52, 860 (2009). 
[21] U. Sivan and Y. Imry, Phys. Rev. B 33, 551 (1986). 
[22] P. N. Butcher, J. Phys.: Condens. Matter 2, 4869 (1990). 



